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Abstract 

The Magnetic Eden Model (MEM) [N. Vandewalle et al., Phys. Rev. E. 
50, R635 (1994)] with ferromagnetic interactions between nearest-neighbor 
spins is studied in {d + 1)— dimensional rectangular geometries for d = 1,2. 
In the MEM, magnetic clusters are grown by adding spins at the boundaries 
of the clusters. The orientation of the added spins depends on both the ener- 
getic interaction with already deposited spins and the temperature, through 
a Boltzmann factor. A numerical Monte Carlo investigation of the MEM has 
been performed and the results of the simulations have been analyzed using 
finite-size scaling arguments. As in the case of the Ising model, the MEM in 
d = 1 is non-critical (only exhibits an ordered phase at T = 0). In d = 2 the 
MEM exhibits an order-disorder transition of second-order at a finite tem- 
perature. Such transition has been characterized in detail and the relevant 
critical exponents have been determined. These exponents are in agreement 
(within error bars) with those of the Ising model in 2 dimensions. Further sim- 
ilarities between both models have been found by evaluating the probability 
distribution of the order parameter, the magnetization and the susceptibility. 
Results obtained by means of extensive computer simulations allow us to put 
forward a conjecture which establishes a nontrivial correspondence between 
the MEM for the irreversible growth of spins and the equilibrium Ising model. 
This conjecture is certainly a theoretical challenge and its confirmation will 
contribute to the development of a framework for the study of irreversible 
growth processes. 



1 Introduction 



The study of kinetic growth models such as directed percolation, Eden growth, bal- 
listic deposition, diffusion limited aggregation, random deposition with and without 
relaxation, cluster-cluster aggregation, etc., is motivated by their interest in many 
areas of scientific research and technology such as polymer science, crystal and poly- 
crystalline growth, gelation, fracture propagation, epidemic spreading, bacterial and 
fungi growth colonies, colloids, etc. [1, 2, 3, 4, 5]. Within this context the Eden 
model [6] has become an archetype growth model. Eden clusters are compact but 
the self-affinity that characterizes the behavior of the growing interface is of much 
interest (see e.g. [7, 8, 9, 10, 11, 12]). Few years ago Ausloos et al. [13] have intro- 
duced an additional degree of freedom to the Eden model, namely the spin of the 
added particles. More recently, the Eden growth of clusters of charged particles has 
also been studied [14]. 

In the Magnetic Eden Model (MEM) [13] with spins having two orientations (up 
and down) the growth of the cluster starts from a single seed, e.g. a spin up seed, 
placed at the center of the two-dimensional square lattice, whose sites are labelled 
by their rectangular coordinates (i,j). Then, the growth process of the resulting 
magnetic cluster consists in adding further spins to the growing cluster taking into 
account the corresponding interaction energies. By analogy to the Ising model [15] 
one takes J as the coupling constant between nearest-neighbor (NN) spins Sij and 
the energy E is then given by 

E S^,S^i'^ (1) 
(U>i'j'> 

where means that the summation is taken over occupied NN sites. The 

spins can assume two values, namely S^j = ±1. Throughout this work we set the 
Boltzmann constant equal to unity (fc^ = 1)) we consider J > (i.e., the ferromag- 
netic case) and we take the absolute temperature T measured in units of J. In the 
MEM a spin is added to the cluster with a probability proportional to the Boltz- 
mann factor exp(— AE/T), where AE is the total energy change involved. It should 
be noted that at each step all sites of perimeter are considered and the probabilities 
of adding up and down spins have to be evaluated. After proper normalization of 
the probabilities the growing site and the orientation of the spin are determined 
through a pseudo-random number generator. 

It is worth mentioning that the MEM has originally been motivated by the study 
of the structural properties of magnetically textured materials [13]. While these pre- 
vious studies of the MEM were mainly devoted to determine the lacunarity exponent 
and the fractal dimension of the set of parallel oriented spins [13], the aim of the 
present work is to complement these previous investigations by studying the critical 
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behavior of the MEM using extensive Monte Carlo simulations and applying a finite- 
size scaling theory. Also, our study is performed in confined (stripped) geometries 
which resemble recent experiments where the growth of quasi-one-dimensional strips 
of Fe on a Cu(lll) vicinal surface [16] and Fe on a W(llO) stepped substratum [17] 
have been performed. In fact, the preparation and characterization of magnetized 
nanowires is of great interest for the development of advanced microelectronic de- 
vices [16, 17, 18, 19]. Furthermore, the growth of metallic multilayers of Ni and Co 
separated by a Cu spacer layer has recently been also studied [20] . 

Another goal of the present work is to compare the results obtained for the MEM 
with the well known behavior of the classical Ising model [15, 21], an archetypical 
model in the study of thermally driven (reversible) phase transitions in equilibrium 
systems. The Ising Hamiltonian (H) is given by 

H = SijS.,, , (2) 

where means that the summation runs over all NN sites, Sij = ±1 is the 

state of the spin at the site of coordinates and J is the coupling constant ( 
J> 0). 

The MEM is also similar to a family of models for the stochastic growth of 
crystals generically known as crystal growth models (CGM) [22, 23, 24, 25, 26], for 
a review see e.g. [27]. As in the MEM, in the case of CGM each atom is adsorbed 
with a given probability conditional to the actual configuration of neighboring atoms 
on the previous layer (s). However, in contrast to the MEM, the crystal is supposed 
to grow layer after layer. It should also be noticed that relationships established 
between CGM and a special class of Ising models [24, 26, 28] have allowed to derive 
exact results. Therefore, useful comparisons with the MEM will be also discussed 
in the presentation of our results. 

This paper is organized as follows: in Section 2 we give details on the simula- 
tion method. Section 3 is devoted to the presentation and discussion of the results 
obtained for the MEM in (1 -|- l)-dimensions, while Section 4 refers to results corre- 
sponding to (2 -|- l)-dimensions. In Sections 3 and 4, detailed discussions comparing 
our results with the behavior of the Ising magnet are outlined. Finally our conclu- 
sions are stated in Section 5. 

2 Description of the simulation method 

The MEM in (1 -|- 1)— dimensions is studied in the square lattice using a rectangular 
geometry L x M with M ^ L and imposing periodic boundary conditions along 
the L— direction. The location of each site on the lattice is specified through its 
rectangular coordinates (1 < i < M, 1 < j < L). The starting seed for the 
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growing cluster is a column of parallel oriented spins placed at i = 1. It should 
be noticed that previous simulations of the MEM were restricted to rather modest 
cluster sizes, i.e. containing up to 8000 spins [13], while in the present work clusters 
having up to 10^ spins have been typically grown. We have also studied the MEM 
in (2 + 1)— dimensions employing a L x L x M geometry (M ^ L) with periodic 
boundary conditions along both L— directions. 

3 Study of the MEM in (1+1)— dimensions: results 
and discussion 

Magnetic Eden clusters grown on a stripped geometry of finite linear dimension 
L at low temperatures show an interesting behavior that we call magnetization 
reversal. In fact, we have observed that long clusters are constituted by a sequence 
of well ordered magnetic domains. Spins belonging to each domain, of average 
length Id ^ L, have mostly the same orientation and consecutive domains have 
opposite orientation. Let Ir be the characteristic length for the occurrence of the 
magnetization reversal. Since Ir ~ L, we then conclude that the problem has two 
characteristic length scales, namely Id and Ir, such that Id ^ Ir- 

In ordinary thermally driven phase transitions, the system changes from a disor- 
dered state at high temperatures to a spontaneously ordered state at temperatures 
below some critical value Tc where a second-order phase transition takes place. Re- 
garding the Ising model, one has that, in the absence of an external magnetic field 
{H — 0), the low temperature ordered phase is a state with non- vanishing sponta- 
neous magnetization [±Msp). This spontaneous symmetry breaking is possible in 
the thermodynamic limit only. In fact, it is found that the magnetization M of a 
finite sample can pass with a finite probability from a value near +Msp to another 
near —Mgp, as well as in the opposite direction. Consequently, the magnetization 
of a finite system, averaged over a sufficiently large observation time, vanishes at 
every positive temperature. The equation M (T, = 0) ~ holds if the observation 
time (tobs) becomes larger than the ergodic time (terg), which is defined as the time 
needed to observe the system passing from iM^p to ^Msp. Increasing the size of 
the sample the ergodic time increases too, such that in the thermodynamic limit 
ergodicity is broken due to the divergence of the ergodic time, yielding broken sym- 
metry. Since Monte Carlo simulations are restricted to finite samples, the standard 
procedure to avoid the problems treated in the foregoing discussion is to consider 
the absolute magnetization as an order parameter [29]. Turning back to the MEM, 
we find that the phenomenon of magnetization reversal also causes the magnetiza- 
tion of the whole cluster to vanish at every non-zero temperature, provided that the 
length of the cluster Iq (which plays the role of tabs) is much larger than Id (which 
plays the role of terg )■ Therefore, we have measured the mean absolute column 
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Figure 1: Plots of the probability distribution of the mean column magnetization 
Pi(m) versus m for the fixed lattice width L = 128 and different temperatures, 
as indicated in the figure. The sharp peaks at m = ±1 for T = 0.45 have been 
truncated, in order to allow a detailed observation of the plots corresponding to 
higher temperatures. This behavior resembles that of the one-dimensional Ising 
model. More details in the text. 

magnetization, given by 

1 ^ 

\m{i,L,T)\ = j\Y.Sij\ . (3) 
^ j=i 

In the stripped geometry used in this work the bias introduced by the lineal 
seed (a starting column made up entirely of up spins) can be avoided by calculating 
relevant properties after disregarding spins within a distance approximately equal 
to few times L from the seed. The procedure of column averaging out from the 
transient region represents a significant advantage of the stripped geometry used for 
the simulation of the MEM. In fact, when a single seed at the center of the sample 
is used, the definition of the average magnetization of the whole cluster is strongly 
biased by the cluster's kernel orientation at the early stages of the growing process. 
In addition, using several randomly generated seeds we could also establish that the 
system evolves into a given stationary state independently of the seed employed. 

The mean column magnetization is a fluctuating quantity that can assume L + 
1 values. Then, for given values of both L and T, the probabihty distribution 
of the mean column magnetization (Pi(m)) can be evaluated, since it represents 
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the normalized histogram of m taken over a sufficiently large number of columns 
in the stationary region [30, 31, 32]. In the thermodynamic limit the probability 
distribution {P^{m)) of the order parameter of an equilibrium system at criticality 
is universal (up to rescaling of the order parameter) and thus it contains very useful 
and interesting information on the universality class of the system [33, 34, 35]. 
Figure 1 shows the thermal dependence of PiiiTi) for a fixed lattice size (L = 128) 
as obtained for the MEM. At high temperatures Pl(*^) is a Gaussian centered 
at m = but when the temperature gets lowered, the distribution broadens and 
develops two peaks at m = 1 and m = —1. Further decreasing the temperature 
causes these peaks to become dominant while the distribution turns distinctly non- 
Gaussian, exhibiting a minimum just at m = 0. The emergence of the maxima at 
m = ±1 is quite abrupt. This behavior reminds us the order parameter probability 
distribution characteristic of the one dimensional Ising model. In fact, for the well 
studied d— dimensional Ising model [32, 36], we know that for T > Tc, Pl{M) is a 
Gaussian centered at M = 0, given by 

P^M) oc exp (^^^ j , (4) 

where the susceptibility x is related to order parameter fluctuations by 

Decreasing temperature the order parameter probability distribution broadens, 
it becomes non-Gaussian, and near Tc it splits into two peaks that get the more sep- 
arated the lower the temperature. For T < Tc and linear dimensions L much larger 
than the correlation length ^ of order parameter fluctuations, one may approximate 
Pl{M) near the peaks by a double-Gaussian distribution, i.e. 

Pl{M) oc exp J + exp J , (6) 

where Msp is the spontaneous magnetization, while the susceptibility x is now given 

by 

X = ^ ((M^) - ■ (7) 

Prom equation (4) it turns out that the Gaussian squared width associated with 
high temperature distributions is very close to the 2nd moment of the order param- 
eter, i.e. 

(7^ f« (M^) . (8) 
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Figure 2: Data for strip widths in the range 16 < L < 1024 as indicated in the 
figures, (a) In-hnear plots of L(m^) versus T~^. The slope of the solid fine (linear 
fit to the data) is a = 1.6. (b) Plots of (|m|) versus T. More details in the text. 
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Equation (8) is a consequence of the Gaussian shape of the order parameter prob- 
abihty distribution and, thus, it holds for the MEM as welL From the known 
one-dimensional exact solution for a chain of L spins [37, 38] one can obtain 

X = ^exp(2/r) (9) 

then, equations (5) and (9) lead us to 

(M2) = iexp(2/r) (10) 

(where it has been taken into account that (M) = due to finite-size effects, ir- 
respective of temperature). Prom equations (8) and (10) we can see that the high 
temperature Gaussian probability distribution broadens exponentially as T gets low- 
ered, until it develops delta-like peaks at M = ±1 as a consequence of a boundary 
effect on the widely extended distribution. It should be noticed that for d > 2 
this phenomenon is prevented by the finite critical temperature which splits the 
Gaussian, as implied by equation (6). 

Turning back to the MEM, figure 1 strongly suggests that an analogous mech- 
anism should be responsible for the thermal dependence exhibited by the MEM's 
order parameter distribution function. So, by analogy to equation (9), we assume 
the relation ^ 

X=^exp(a/T) (11) 

to hold for the MEM, where we have introduced a phenomenological parameter a, 
and the susceptibility x is given by equation (5). We find an excellent agreement 
to the data by choosing the value a — 1.6 as observed in figure 2(a), where In- 
linear plots of L{m'^) versus 1/T are shown for strip widths varying in the range 
16 < L < 1024. Figure 2(b) shows plots of versus T for the same lattices. 

This figure shows that increasing L the order parameter curves approach the one 
that corresponds to the thermodynamic limit (i.e., (|m|) = 0{T), where 9 is the 
Heaviside function). 

However, it should be pointed out that the results obtained for the (1 + 1)-MEM 
and the 1-dimensional Ising model do not exactly coincide for finite lattices, as figure 
3(a) shows for the case of the magnetization. Anyway, this fact should not alarm 
us, since it can be seen that differences in the results obtained for both models are 
a direct consequence of the finite-size nature of the lattices used in the simulations 
and consequently they tend to vanish in the thermodynamic limit. This is actually 
shown by figure 3(b), where log-linear plots of {\M\fsing){L,T) — {\m\MEM){L,T) 
versus for two different fixed values of temperature are presented. Thus, we 
conclude that in view of the full qualitative and quantitative agreement between 
both models we can safely estabhsh that, as in the Id Ising model, the (1 -|- 1)-MEM 
is not critical (i.e. it also undergoes a phase transition at Tc — 0). 
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Figure 3: Comparison of results corresponding to the (1 + 1)— MEM and the d = 1 
Ising model, (a) Plots of (|m|) versus T obtained for a lattice of side L — 128. 
(b) Linear-log plots of {\M\i){L,T) - {\m\MEM){L,T) versus L''^ for T = 0.5 and 
T — 1.0. Hence, differences in the magnetization due to finite-size effects appear to 
vanish in the thermodynamic limit. 
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Figure 4: Plots of P{nNN) versus T for unn = 1,2,3, 4 as indicated in the figure. 
The hues are guides to the eye. More details in the text. 

We have also computed the number of already occupied NN sites every time 
that a new particle was added to the spin system, and thus we have obtained the 
normalized probability P{nNN) of having unn occupied NN sites. Figure 4 shows 
the behavior of P{nNN) as a function of temperature. Using this probability we 
have evaluated {unn) — 2.0000(1) irrespective of the temperature. This result can 
be understood considering that the growing process that leads to the assignment of 
a spin Sij = ±1 to each lattice site of coordinates can be studied by means 

of a bond model. In fact, we can assign a bond to each pair of neighboring sites, 
pointing from the earlier occupied site to the later occupied one. So, the process 
that leads to a given spin configuration can be specified by the fields bu{i,j) and 
bR{i,j), where sub-indexes U and R refer to the upper bond of (i, j) (i.e., the bond 
that connects the site of coordinates with that of coordinates + 1)), and to 
the right bond of the site of coordinates respectively [39]. We take b{i,j) = +1 
if the bond points outwards and b{i,j) = —1 if it is directed inwards. Therefore, the 
net bond flux at a given lattice site is given by: 

(f>{i,j) - bu{i,j) + bniij) - bu{i,j - 1) - 6r(z - 1, j) (12) 

and the possible values that can take are = —4, —2, 0, 2. After some algebra, 
it follows that = |(4 — 0) holds for every site on the lattice. Moreover, it 
can be seen that, for an arbitrary d-dimensional lattice of coordination number q, 
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Figure 5: Plots of P{AE) versus T for AE = 0, ±1, ±2, ±3, ±4 (in units of J) as 
indicated in the figure. The hues are guides to the eye. More details in the text. 

nNN = |(g - 0). Then, 

(riArjv) = 2 ^ 

is the mean number of occupied NN, since (0) = 0. For the two-dimensional square 
lattice, q — A and equation (13) yields {unn) = 2, in agreement with the result we 
have already obtained by means of Monte Carlo simulations. 

Further insights into the MEM's growing process can be gained by studying the 
mean energy change involved in the addition of a new particle to the system. The 
process of adding a new spin involves an energy change AE and from the definition 
of the (1 + 1)-MEM, the possible values that AE can take are 0, ±1, ±2, ±3, ±4 
(in units of J). Figure 5 shows plots of the normalized probabihty P{AE) versus 
T for each of these values. The non-equilibrium natiirc of the MEM manifests itself 
through much more complex probability distributions P{AE) (see figure 5) than the 
corresponding to the equilibrium Ising model where AE can take only five different 
values, namely 0, ±2, ±4 (in units of J). The results shown in figures 4 and 5 confirm 
the non-trivial nature of the fink established between the MEM at stationarity and 
the Ising model in equilibrium. 

It should be noticed that for the case of crystal growth models (CGM) [22, 23, 24] 
the growing conditions are quite different than those of the MEM. In fact, in CGM 
the crystal grows layer by layer in a given direction [24, 26]. Furthermore, the prob- 
abihty distribution of the predecessor spin layer is sampled from the equilibrium 
distribution, so will be the probability of spins in subsequent layers. This particular 
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Figure 6: Plots of the probability distribution PL(m') versus m for the fixed lattice 
size L = 16 and different temperatures, as indicated in the figure. The occurrence of 
two maxima located at m' — iM^p (for a given value of Mgp such that < Mgp < 1) 
is the hallmark of a thermal continuous phase transition that takes place at a finite 
critical temperature. 

growth mechanism allow to establish dual transformations with the kinetic Ising 
model [25, 26] and to extract some exact results. In contrast, the growing inter- 
face of the MEM is self-affine and the system is far from equilibrium. So, the link 
between the Id Ising model and the (1 + 1)— dimensional MEM is quite challenging. 

4 Study of the MEM in (2+1)— dimensions: results 
and discussion 

4.1 The order parameter and its probability distribution 
function 

In order to compare the (2 + 1)— dimensional MEM and the 2d Ising model, we have 
first studied the order parameter probability distribution PL{m ), where m takes 
now + 1 possible values (see figure 6). For high temperatures, the probabihty 
distribution corresponds to a Gaussian centered at m =0. At lower temperatures 
we observe the onset of two maxima located at m = itM^p (0 < Mgp < 1), which 
become sharper and approach m = ±1 as T is gradually decreased. 
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Figure 7: Plots showing the location of the maximum of the probability distribution 
as a function of temperature for both (d+ 1)— dimensional MEM models {d — 1,2). 
The lines are guides to the eye. The smooth transition for d=2 constitutes another 
evidence of the finite critical point associated with the (2 + 1)— MEM. 



Figure 7 shows the location of the maximum of the probability distribution as a 
function of temperature for both [d + l)-dimensional MEM models, with d = 1,2. 
While for the d — 2 case we observe a smooth transition from the m^^^ — value 
characteristic of high temperatures to nonzero m^^^ values that correspond to lower 
temperatures, the curve obtained for = 1 shows, in contrast, a Heaviside-like jump. 
As already discussed, the behavior exhibited by the (2 + l)-dimensional MEM (e.g., 
as displayed by figures 6 and 7) is the signature of a thermal continuous phase 
transition that takes place at a finite critical temperature. 

The broken symmetry at a finite critical temperature implied by the thermal 
continuous phase transition can be explained in terms of the broken ergodicity that 
occurs in the system when we tend to the thermodynamic limit (L oo) making 
use of the temperature dependence exhibited by the order parameter distribution 
function. In fact, if we set the characteristic length of MEM's domains Id equal to 
an ergodic length lerg, we can carry out a complete analogy with the Ising model by 
associating l^rg to t^rg (the Ising model ergodic time) and the cluster's total length 
Ic (already defined in Section 3) to the Ising model observation time tof,^. In this 
way, we encounter that excursions of m from m = +Msp to m = —Mgp and vice 
versa occur at length scales of the order of lerg- When the cluster's total length 
becomes larger and larger [Ic ^ Urg) the whole cluster's magnetization is averaged 
to zero. Furthermore, l^rg diverges as the strip's width becomes larger and larger. 
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and again broken symmetry arises as the consequence of broken ergodicity. 

4.2 Order-disorder phase transition in the (2+l)-dimensional 
MEM: Finite-size effects and scahng analysis 

As already anticipated and as it follows from figures 6 and 7, the (2 + l)-dimensional 
MEM exhibits a thermally driven order- disorder transition at a finite temperature. 
In the thermodynamic limit (L — > oo) we expect to determine a critical temperature 
Te such that (|m |) = for T > while (|m |) remain non- vanishing at temperatures 
below Te. 

From the finite-size scaling theory, developed for the treatment of finite-size ef- 
fects at criticality and under equilibrium conditions [40, 41], it is well known that 
if a thermally driven phase transition occurs at a temperature Tc > in the ther- 
modynamic limit, then in a confined geometry this transition becomes smeared out 
over the temperature region AT(L) around a shifted effective transition temperature 
Tc{L), and the following relationships hold: 

AT{L) oc L-^ , (14) 

and 

|T,(L) - Tel oc L-^ , (15) 

where the rounding and shift exponents arc given by ^ = A = z/~^, respectively 
(recalling that u is the exponent that characterizes the divergence of the correlation 
length at criticahty). 

Furthermore, from well established finite-size scahng relations, the following 
Ansatze hold just at criticahty {T — Tc) : 

(|m'(L,r = Te)|)ocL-^/'^ (16) 

and 

Xma.m^V''' , (17) 

where j3 and 7 are the order parameter and the susceptibility critical exponents, 
respectively. Note that Xmax{L), as given by equation (17), refers to the maximum 
of x(L, T) as a function of T for fixed lattice size L. 

In view of the encountered analogies between the MEM and the Ising model, 
it is natural to test the validity of equations (14-17) for the case of the MEM in 
(2 + l)-dimensions. It should be noted that as in the case of equilibrium systems, 
in the present case various "effective" L-dependcnt critical temperatures can also 
be defined. In particular, we will define Tci(L) as the value that corresponds to 
(|m I) = 0.5 for fixed L, and Tc2{L) as the one corresponding to the maximum of 
the susceptibility for a given L, assuming that the susceptibility is related to order 
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Figure 8: Plots of Tcn{L) versus L ^ (for n = 1, 2). The solid lines show the linear 
extrapolations that meet at the critical point given by Tc — 0.69 ± 0.01. 

parameter fluctuations in the same manner as for equilibrium systems (as given by 
equations (5) and (7)). Then, we should be able to obtain Tc from plots of Tc„(L) 
versus (for n = 1, 2), as it is shown in figure 8. Following this procedure we find 
that both Tci(L) and Tc2{L) extrapolate (approximately) to the same value, allowing 
us to evaluate the critical temperature T^. — 0.69 ±0.01 in the thermodynamic limit. 

After determining T^, the correlation length exponent u can be evaluated by 
means of equation (15), making the replacement \ = l/u. In fact, taking at 
the mean, maximum and minimum values allowed by the error bars, wc obtain six 
log- log plots of \Tcn{L) — Tc\ versus L for n = 1, 2. The slope of each of these plots, 
not shown here for the sake of space, yields a value for u. The obtained values are: 

i/ = 1.08(7; = 0.68), u = 1.00(rc = 0.69), 

i/ = 0.88(rc = 0.70) for n = 1, (18) 

and 

u = 1.20(7; = 0.68), u = 1.08(7; = 0.69), 

1/ = 0.95 (Tc = 0.70) for n = 2. (19) 

Thus our estimate is given hy = 1.04 ± 0.16, where the error bars refiect the 
error derived from the evaluation of Tc, as well as the statistical error. 

Figure 9 shows plots of the susceptibility versus T as obtained using lattices of 
different side. It is found that the susceptibility exhibits a peak which becomes 
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Figure 9: Behavior of the susceptibihty as a function of temperature. Each curve 
shows a peak which becomes sharper and shifts towards lower temperatures as L is 
increased. The inset shows in greater detail the peaks corresponding to the smaller 
lattices (L = 16,24,32). 

sharper and shifts towards lower temperatures when L is increased. This behavior 
is, in fact, already anticipated by equation (17), and it allows us to evaluate 'j/u 
from the slope of a log-log plot of Xmax versus L, as figure 10 shows. The linear fit 
yields ^/v — 2.02 ± 0.04. Using this value and the value formerly obtained for v we 
thus determine 7 = 2.10 ± 0.36. 

Figure 11 shows log-log plots of (|m |)(T = T^) versus L for the mean, maximum 
and minimum allowed values of T^. Considering only the larger lattices, the linear 
fits to the data according to equation (16) yield the following estimates: fijv — 0.11, 
i5jv — 0.16 and — 0.19. We then assume the value = 0.15 ± 0.04, where 
the error bars reflect the error derived from the evaluation of Tg, as well as the 
statistical error. From this value and the value formerly obtained for v we thus 
determine /3 = 0.16 ± 0.05. 

The critical exponents of the MEM in (2 -|- 1)— dimensions as obtained using 
a finite-size scaling analysis are, so far: v — 1.04 ± 0.16, 7 = 2.10 ± 0.36 and 
/9 = 0.16 ± 0.05. If we recall the exactly known critical exponents of the d = 2 Ising 
model, i.e. z/ = 1, 7 = 7/4 and /5 = 1/8, we find that the (2 + 1)— MEM has the same 
critical exponents within error bars. These results further support our conjecture 
on the connection between the MEM in (2 -|- l)-dimensions and the Ising model in 
2 dimensions. 

As in the case of the MEM ind — 1, we have also computed the number of already 
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Figure 10: Log-log plot of Xmax versus L. The linear fit (solid line) yields ^/v — 
2.02 ± 0.04. 
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Figure 11: Log-log plots of {m {T = Tc)) versus L ^ for the mean, maximum and 
minimum allowed values of Tc. The linear fits (solid lines) yield an estimate [3/u — 
0.16 ±0.05. 
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occupied NN sites every time that a new particle was added to the spin system. We 
found that the value (unn) = 3.0000(1) holds for all temperatures, which is indeed 
the result given by equation (13), since g = 6 for the three-dimensional square 
lattice. 

At this stage, we may recall that for the (1 + 1)— MEM (njvTv) equals the co- 
ordination number of the d = 1 Ising model, and that we found that both models 
have the same critical temperature and exhibit the same critical behavior. Rea- 
soning by analogy, we may expect a coincidence between the critical temperature 
for the (2 -|- 1)— MEM and the corresponding one for a = 2 Ising model defined 
on a lattice of coordination number g = 3. However, this comparison cannot be 
carried out, since the critical temperature of an Ising model depends on both the 
coordination number q and the topological structure of the lattice, but for o? > 2 and 
a given value of q the topological structure is not unique. For instance, for d 2 
and q = 3, we can pass from the honeycomb lattice (HL) to the expanded Kagome 
lattice (EKL) through the application of a star-triangle transformation and obtain 
the exact values of their critical points, which turn out to be [42] Tc = 1.5187 (HL) 
and Tc = 1.4530 (EKL). 

5 Conclusions 

In the present work we have studied the growth of magnetic Eden clusters with 
ferromagnetic interactions between nearest neighbor spins in a (d+ 1)— dimensional 
rectangular geometry (for d = 1,2), using Monte Carlo simulations and applying a 
finite-size scaling theory. The results obtained allow us to conjecture a nontrivial 
correspondence between the MEM for the irreversible growth of magnetic materials 
and the classical Ising model under equilibrium conditions. In fact, we have found 
that the {d + 1)— dimensional MEM and the (i— dimensional Ising model behave 
identically (unless finite-size differences that vanish in the thermodynamic limit) at 
criticality, i.e. that both models belong to the same universality class. We also con- 
jecture that this correspondence would remain at higher dimensions (li > 2). The 
results obtained strongly suggest a link between the temporal evolution of equilib- 
rium systems and the stationary growth of nonequilibrium systems. We thus believe 
that this work will stimulate further developments in the field of nonequilibrium ki- 
netic growth models. A more precise numerical test of the posed conjecture will 
certainly require a considerable computational effort but it will be of great interest. 
Furthermore, analytical developments aimed to establish a theoretical framework 
for the understanding of far from equilibrium growth phenomena will become stim- 
ulated by the reported findings. 
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